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ON THE BACKWARD BEHAVIOR OF SOME DISSIPATIVE 
EVOLUTION EQUATIONS 


YANQIU GUO AND EDRISS S. TITI 


Abstract. We prove that every solution of a KdV-Burgers-Sivashinsky type equa¬ 
tion blows up in the energy space, backward in time, provided the solution does not 
belong to the global attractor. This is a phenomenon contrast to the backward 
behavior of the periodic 2D Navier-Stokes equations studied by Constantin-Foias- 
Kukavica-Majda [T8], but analogous to the backward behavior of the Kuramoto- 
Sivashinsky equation discovered by Kukavica-Malcok [50]. Also we study the back¬ 
ward behavior of solutions to the damped driven nonlinear Schrodinger equation, 
the complex Ginzburg-Landau equation, and the hyperviscous Navier-Stokes equa¬ 
tions. In addition, we provide some physical interpretation of various backward 
behaviors of several perturbations of the KdV equation by studying explicit cnoidal 
wave solutions. Furthermore, we discuss the connection between the backward be¬ 
havior and the energy spectra of the solutions. The study of backward behavior 
of dissipative evolution equations is motivated by the investigation of the Bardos- 
Tartar conjecture stated in j5]. 


1. Introduction 


Consider a KdV-Burgers-Sivashinsky (KBS) type equation on the torus T = [—§,■§]: 


( 1 . 1 ) 


Ut T ttUx (3u -(“ /*, t G T, 


with the initial condition w(0) = Uq. We assume that v > 0, /3 G R, and 7 ^ 0. Also, 
we suppose that the time independent forcing / and the initial data u 0 both have 
spatial mean value zero, i.e., f T fdx = f T u 0 dx = 0, and hence f T u(x, t)dx = 0 for all 
time t within the lifespan of the solution u. 

We define the spaces 
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We call a solution u(t) of (II.ID global if u(t) is defined for all t 6 1. The main 
result of this paper states that, if v > 0, then there does not exist a global solution 
of (11.11) in the energy space H unless u belongs to the global attractor. Specifically, 
we have the following theorem. 

Theorem 1.1. Consider the following two cases: 

(i) Assume v = 0 and f G H . For any uq G H, \1.1\) has a unique global solution 
u(t ) G H for all t G I. 

(ii) Assume v > 0 and f G V'. Let u(t) : [0, oo) -A H be a solution of M.l\) which 
does not belong to the global attractor. Then u{t) cannot be extended to a global 
solution for all t G R. 

The study of the behavior of solutions to (11.11) backward in time is motivated by 
the pioneering work [IB], where Constantin, Foias, Kukavica and Majda investigated 
the backward behavior of the solution to the 2D periodic Navier-Stokes equations 
(NSE). Indeed, they showed that the set of initial data for which the solution exists 
for all negative time and has exponential growth is rather rich (dense in FL with the 
topology of (Hp er )' where FL = {u G L^ er { T 2 ) 2 : divw = 0 ,f T2 udx = 0}), which is 
a quite remarkable result since it indicates that the backward behavior of the 2D 
periodic NSE is closer to the corresponding linear dissipative equation. Also their 
result provided a partial positive answer to a conjecture in [5], where Bardos and 
Tartar conjectured that, for the 2D periodic NSE, the solution semi-flow S[t)FL is 
dense in the phase space FI for any fixed t > 0. On the other hand, the 2D Euler 
equations are globally well-posed forward and backward in time, and the energy is 
conserved for smooth solutions. The result in [T8j may be rephrased as, if we add 
a viscosity to the 2D periodic Euler equations, which gives the NSE, then there is 
still a rich set of initial data such that the solution can be extended to all negative 
times with the energy growing exponentially, i.e., adding a viscosity to the periodic 
2D Euler equations fails to drive these solutions to blow up backward in finite time. 

We remark that a corollary of Theorem ll.il (i.e. Corollary 12.71 in section [2]) implies 
that the solution semi-flow S(t)H of the KBS equation (11.11) with v > 0 is not dense 
in FI for any t > 0. This result is contrast to the Bardos-Tartar conjecture on the 
NSE mentioned above [5]. 

Notice that, if we remove the viscosity from the KBS equation (II. ip . i.e. v = 0, 
then it becomes 

u t uu x /3u + r yu xxx = f , iGT, (1.2) 

which is a KdV type equation with a linear term — (3u. It is well-known that the 
real-valued KdV equation is globally well-posed in H for t G R (see, e.g., mm, and 
references therein; see also [38] for a similar result concerning a coupled system of 
KdV). If f3 < 0, then forward in time, (11.21) is a weakly damped KdV equation that 
possesses a global attractor in H (33j[36j. On the other hand, if [3 > 0, then — /3u plays 
the role of a source term forward in time as well as a weak damping backward in time. 
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If we consider (1 1.2 j) backward in time with f3 > 0, it is a dissipative system possessing 
a global attractor, which seems more “stable” than 2D periodic Euler equations. So 
it could be reasonable to guess that the solution of the KBS equation (II. ip (i.e. (\\.2\ 
plus a viscosity) might have a similar behavior backward in time as the 2D periodic 
NSE. However, such intuition turns out to be completely false. Indeed, Theorem 
11.11 tells us that no solutions of (11. ip . with v > 0, can be extended to all negative 
time unless they belong to the global attractor. Here, we would like to remark that 
the different backward behavior of the KBS equation (11.11) and the 2D periodic NSE 
might be related to the rate of the energy cascade to small scales for the Burgers and 
the KdV equations as well as for the Euler equations, respectively. Furthermore, we 
stress that, for the backward blow-up of the KBS equation (II. ip . the viscosity plays a 
vital role, contrast to the backward blow-up of the viscous Burgers equation. Indeed, 
the inviscid Burgers equation can also blow up backward in finite time, but the KBS 
equation without viscosity (i.e. the KdV type equation (II. 2p ) is globally well-posed 
for all t G M. 

The proof of Theorem 11.11 follows the idea presented in [50]. It is shown in [50] 
that all solutions of the Kuramoto-Sivashinsky equation off the global attractor blow 
up backward in finite time. Notice that, the Kuramoto-Sivashinsky equation involves 
a hyperviscosity tending to drive the energy increasing very fast backward in time. 
However, the KBS equation (II.ip contains a regular viscosity, and a dispersion dis¬ 
tributing the energy, which seems closer to the NSE. Nonetheless, the energy of the 
KBS equation still increases too fast and blows up backward in finite time, like the 
Kuramoto-Sivashinsky equation, due to the negative viscosity. It can be viewed as 
adding the dispersion term u xxx to the viscous Burgers equation fails to prevent so¬ 
lutions blowing up as the time goes backwards. 

In addition to the study of the KBS equation (II.ip . we investigate the backward be¬ 
havior of solutions to some other typical dissipative equations. Specifically, we show 
that any solution off the global attractor, to the damped driven nonlinear Schrodinger 
equation, grows exponentially fast backward in time. Also, we prove that all solu¬ 
tions to the complex Ginzburg-Landau equation of certain parameter regime, blow 
up backward in finite time, unless they belong to the global attractor. Moreover, 
we discuss the asymptotic behavior of solutions to a hyperviscous NSE, forward and 
backward in time. 

The viscous KdV equation is a special case of the KBS equation by taking f3 = 
7 = 0 and / = 0 in (II. ip . It is well-known that the classical periodic KdV equation 
has explicit cnoidal wave solutions. Also the solutions of the KdV are connected 
to the spectrum problem of certain Schrodinger operator. These nice properties of 
the KdV equation provide us an opportunity to employ the modulation asymptotic 
analysis to study the backward behavior of the viscous KdV equation. Indeed, we 
give an evidence that the amplitude of the cnoidal wave solution of the viscous KdV 
equation approaches infinity backward in finite time, which is consistent with our 
result on the backward blow-up of the energy. Furthermore, we attempt to analyze 
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the underlying mechanism for the backward blow-up of the viscous KdV equation 
as well as the backward non-blowup phenomenon of the 2D periodic NSE, by using 
the energy spectra of solutions together with the Kolmogorov theory of turbulence. 
These discussions from the physical point of view are novelties of the manuscript and 
are major motivation for writing this article. 

Now we briefly survey some related works in the literature concerning the backward 
behaviors of dissipative evolution equations. It is proved in [65] that the 2D periodic 
viscous Camassa-Holm equations (also known as Navier-Stokes-ct model [25]) has an 
analogous backward behavior as the 2D periodic NSE. Some different phenomenon 
has been discovered in |T9] for Burgers’ original model for turbulence which consists 
of a PDE coupled with a nonlocal ODE, and it was shown for this system, there are 
three possible behaviors of a solution as t —* — oo: it can be globally bounded, grow 
exponentially fast, or grow faster than any exponential. 

Recall that the 2D periodic NSE reads 

iu t - uAu + (u ■ V)w + Vp = /, xeT 2 

jdivu = 0, 1 ’ 

with u(0) = Uq. A main ingredient in studying the backward behavior of solutions 
to (11.311 is a pair of orthogonality relations of the form ((u ■ X7)v,v ) L 2 = 0 and ((w • 
V)w, Am ) L 2 = 0, provided u, v are periodic in two dimensional space and divergence 
free. Such pair of orthogonality relations lead to a pair of energy and enstrophy 
formula: \f t \u\ 2 L2 + v\ Vu|| 2 = (/, u ) L 2 and \ j- t \ Vu|| 2 + v\ Am|| 2 = —(/, Au) L 2 , which 
are identical to the corresponding estimates for the linear Stokes equations (i.e. the 
equation (II.3ft with the nonlinear term vanished), and are the essence for proving that 
there is a rich set of initial data such that the solution to (II. 3ft exists globally for all 
t 6 K (see [T8] ). It is demonstrated in [26] how to generate a second orthogonality 
property for Lorenz equations, by finding an alternate linear operator. 

In addition, see [20j [66j for some related results on this topic. 

Throughout, | • | and (•,•) represent the norm and the inner product in L 2 er (T) 
respectively, and || • || stands for the norm in H^ er { T), i.e., ||w|| 2 = \u x \ 2 + \u\ 2 for u G 
Hp er (T). As usual, we set A = — d xx . Moreover, we will denote by c a dimensionless 
constant that may change from line to line. 

The paper is organized as follows. In section 11.11 we shall prove the main rigor¬ 
ous mathematical result of the paper: Theorem 11.11 I 11 section [3] we investigate the 
backward behavior of solutions to damped driven nonlinear Schroclinger equation, 
complex Ginzburg-Landau equation, and hyperviscous Navier-Stokes equations, as 
well as mentioning some open problems. Section [4] provides an elegant physical inter¬ 
pretation of the backward behavior of various perturbations of the KdV equation, and 
also includes some interesting discussion about the relation of the backward behavior 
of dissipative equations with the Kolmogorov theory of energy spectra of turbulence. 
The appendix is devoted to proving two results presented in section [3] 
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2. Proof of Theorem 11.11 

To begin with, let us consider the first part of Theorem ll.il i.e., v — 0. Without the 
viscous term, (II. ip becomes a KdV type equation H 1.2 )) . The global wcll-posedness 
of the real-valued KdV equation u t + uu x + u xxx = 0 in H has been established in 
Bourgain’s seminal paper [10] (see ||3] for an alternative proof). Note that the energy 
of the KdV equation is conserved. Adding a linear perturbation —f3u and a time- 
independent forcing / G H to the KdV equation (ending up with equation (ll.2p ) 
destroys the energy conservation, but the global well-posedness in H for all t e I is 
still valid for (II.ip . Indeed, taking the scalar product of (II.2p with u yields 

1 =(/,«). 

By assuming (3 > 0 and using Holder’s, Young’s and Gronwall inequalities, we deduce 
\u(t)\ 2 < e m \u 0 \ 2 + ^(e 3/3t - 1), for t > 0, (2.1) 

and 

|u(t)| 2 < e«|«o| 2 - jy i - 1)|/| 2 , for t < 0, (2.2) 

Notice that, (12.11) implies that the energy grows at most exponentially fast as t —> Too, 
while (12.2p shows that as the time t goes backwards, the energy is uniformly bounded. 
In fact, lim sup 4 _^_ 00 u{t ) < -b |/| 2 , and thus all trajectories of equation (II.2p approach 
a global attractor in H as t —> —oo, provided /3 > 0. By switching the direction of 
the time, the same analysis holds true for the case f3 < 0. 

Now we prove the second part of Theorem 11.11 i.e., the case v > 0. 

2.1. Asymptotic estimates for the long-time dynamics. In the second part of 
Theorem ll.il we assume / G V'. In the proof, for the sake of clarification, we assume 
f3 > 0, which is the more interesting case. The case /3 < 0 can be treated analogously 
with less effort. 

In order to study the behavior of solutions for evolutionary equations forward or 
backward in time, a natural method is to perform an a priori estimate of the energy, 
and for equation (II.ip the energy is the L 2 -norm of the solution. A straightforward 
way of estimating the energy |u| is to take the scalar product of the KBS equation 
(11.11) with u\ 

“T'I m | 2 + u \ u x \ 2 — (3\u\ 2 = (/, u) < — \u x \ 2 + —\A~^f\ 2 . 

2 dt 11 1 u ’ ; - 2 1 1 2v Jl 

It follows that 

d\u\ 2 + (['A, - 2/3)|u| 2 < l|/r4/| 2 , 


( 2 . 3 ) 
















6 


Y. GUO AND E. S. TITI 


where the Poincare inequality \u x \ 2 > Ai|m | 2 has been used, for any u £ V. Here, 
Ai is the first eigenvalue of the operator A = — d xx , i.e., X\ = ( 2-k/L ) 2 . Unless 
the viscosity v is large (i.e., v > which is not an interesting case), the energy 
estimate (12.311 does not provide an uniform bound on the L 2 -norm of u(t) for all 
time t. However, intuitively the energy of the KBS equation (II.ip must be uniform 
bounded since it involves the viscous term vu xx acting as a strong dissipation. Notice 
that, a drawback of the above energy estimate is that it does not take advantage of the 
nonlinear convection term uu x , which plays the role of a medium transferring energy 
from lower to higher wave numbers to prevent the growth of low wave number modes. 
Indeed, in the KBS equation (II.ip . — /3u amplify the lower wave numbers while the 
diffusion vu xx damps the higher wave numbers, and the interaction between these 
two linear terms are through the energy conducting uu x . 

In order to take advantage of the convection uu X) we employ the method of Lya¬ 
punov functions (also called the background flow method (13, [16:21], [22]). More pre¬ 
cisely, we subtract u by an appropriate gauge function 0 and study \u — 0|, instead of 
directly estimating the energy |w|. This technique was used in [Ill i~l. [371 H 8 ] 1501 l58 j 
for the purpose of studying the Kuramoto-Sivashinsky equation (KSE). It is worth 
mentioning that, in the case of odd functions, which is an invariant space for the KSE 
(also for equation ( 11.111 with 7 = 0 ), the gauge function 0 is simpler and was intro¬ 
duced in [5BJ. On the other hand, the best asymptotic energy estimate with respect 
to the length L for the KSE was obtained in [61] (see also [32]) by using different 
techniques. 

We adopt the gauge function 0 which was used in [37] EQ]. Similar idea can be 
found in on Hi- To construct the function 0, we need to do some preparation. 
In fact, for every e G (0,L/2), there exists a periodic non-negative smooth function 
b e : R —* [0, 00 ) with period L and supported in (—e, e) such that 

( 1 ) f T b e (x)dx = L ; 

(2) sup xm b £ (x) < cL/ e; 

(3) | 6 e | < cL/eU 

(4) | 6 '| < cL/et. 

An example of such function was given in [50]: if define rj(x) = L(e — \x\)/e 2 for 
|x| < e, then b e can be obtained by mollifying and periodically extending r/ to the 
whole real line. The following Poincare type inequality was proved in [37] . 

Proposition 2.1. [37] Let b e : R —» [0, 00 ) be a periodic function with period L 
such that suppA fl [— -§,§] C (—e,e), J T b € (x)dx = L, and sup xeIR b e (x) < cL/e. If 
f T b e (x)u(x)dx = 0, for some u e V, then 


b e (x)u 2 (x)dx < cqcL / u 2 x {x)dx, for some cq > 0 . 
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As in [371 HO], we define a periodic smooth function : R —* M by 


rx 

<f>a,e( x ) =oix-a / b e (y)dy, 

Jo 

(2.4) 

where a > 0. Since b e is periodic with period L and f T b e (x)dx = L, we see that 
is also periodic with period L. Furthermore, by using the properties of the function 
b e , simple calculation gives 

\fa,e\ < caLF 

(2.5) 

K,el < ca(L^ + Ljeh), 

(2,6) 

\f"\ < caL/e^. 

(2.7) 


Also the following simple fact is useful: since 4> ai€ is periodic, its L 2 -norm has trans¬ 
lation invariance \<f a ^(x)\ = \4> a ,e(x + 01 f° r an y £ G 

Theorem 2.2. For sufficiently large time t>t\ which depends on |uo|, the solution 
S(t)u 0 of M.l\) enters a ball in H with the radius p satisfying the following asymptotic 
relation with the length L and the parameters v, fd, 7; 

p ~ L 5 / 2 , p ~ —, p ~ fd 5 ^ 2 , p ~ 7. (2.8) 

Remark 2.3. Theorem 12.21 shows that there exists an absorbing ball for (11.11) . so by 
the classical theory of attractors (see, e.g. H). the KBS equation (11.11) possesses a 
global attractor. 

Proof. The approach follows the idea in [37, 50j- We aim to estimate 

F<x,e(t) := inf / (u(x,t) - + f)) 2 dx. 

J t 

To this end, for ( 6 I, we calculate 

“ J (' U{X , t) - f a ,e( X + 0Y dx = ~ f Ut ( X ’ WaA 1 + 0 dx - ( 2 ‘ 9 ) 

The first term on the right-hand side of (12.91) can be estimated by usual way. In fact, 
taking the scalar product of (II.ip with u yields 

— M 2 = -v\u x \ 2 + fd\u\ 2 + (f,u) 

< ~v\u x \ 2 + (d\u\ 2 + i|/H/| 2 , 
by using Cauchy-Schwarz and Young’s inequalities. 


( 2 . 10 ) 
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To calculate the last term on the right-hand side of (12.91) . we take the scalar product 
of (11.11) with 0 Qi<E (x + 0; an d notice that (j>' ae {x) = a — ab e (x), 



< t^KI 2 - (| - P)W\ 2 + A 2/| 2 + cua 2 (L + L 2 / e) 



where we have used the estimate (I2.5|) - (l2.7j) in the last inequality. 


Now, substituting (12 TOR and (12.1 lj) into (12.91) implies 



^ - 2/3)|u| 2 + -\A~if\ 2 + cm 2 (L + L 2 / e ) 
2 v 




( 2 . 12 ) 


Next we intend to apply Proposition 12. II to estimate the last term in (I2.12p . Recall 
that we have defined F a>e (t) = inf^ G T f T ( u ( x ,t) ~ 4>a,e{x + £)) 2 da;. Notice, for each 
£x t, the mapping £ ha f T (u(x,t ) — (p ate (x + £,)) 2 dx is continuous and periodic on K. 
Thus, we can assume the minimum of f T (u(x,t) — (f> a ,e(x + Q) 2 dx occurs at a point 
£*(f) G T. Then £ *(t ) must satisfy 



Recall that is a periodic function, thus J T cp ate (x + ^*(t))(p' a t (x + ^*{t))dx = 0. It 
follows that 


and by (12.41) . this is equivalent to 
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Since we assume u has spatial mean value zero, one has 

j u(x,t)b e (x + £*(t))dx = 0. 

J T 

Consequently, by Proposition 12.II and keeping in mind that u and b e are both periodic 
functions with the period L, we conclude that 


u 2 (x , t)b e (x + t/*(t))dx < cc,eL\u x \ 2 . 


(2.13) 


Now we compute 
/ F v (t) = limsup ~ 


= lim sup 

s—H+ 

< limsup 

s—H+ 

d r 

dt Jt 


f T (u(x, s ) - (p a ,e(x + £*(s))) 2 dx - f T (u(x, t ) - </> a , e (x + C(t))) 2 dx 


s — t 


f T (u(x, s) - 4>a ie (x + £,*(t))) 2 dx - f T (u(x, t ) - 0 a>e (x + £*(t))) 2 oh 


s — t 


(u(x,t) - 4>a,e(x + Ofdx 


£=£*(*) 


< (—— + ac 0 eL)\u x \ 2 — (a — 4:/3)\u\ 2 H —\A ? f\ 2 + cua 2 (L + L 2 /e) 


v 


+ c/3a 2 L 3 + ( —t—a 2 L 2 , 


z/e J 


where (12.12ft and (12.13H were used to obtain the last inequality. 
By choosing a > 8/3, and setting e = . 7a u CoL , we obtain 


d 


Oi, ,2 4 . 1 .2 


t: Fa,e(.t) < — — |u| H —|A 2/| -|- C 0 (a) 


dt 


where C 0 (a) = cva 2 L + ca 3 L 3 + c/3a 2 L 3 + ^a 5 L 5 . 

We shall look for an estimate on F a ^(t) by using (12.141) . Indeed, 


(2.14) 


F afi (t) = inf / (u(x,t) - 4> a ,e{x + i)) 2 dx < 2|u(t)| 2 + 2|</> Qj£ 
£ eT J T 


It follows that 


a 


a 


a 


a 


u{t)\ 2 < ~F ae (t) + ^\4>a,e\ 2 < ~F a ,e(t) + Ca 3 L 3 , 


where (12.51) was used. Substituting the above inequality into (I2.14j) yields 
4 Fa,e(t) < -JF a , e (t) + ca 3 L 3 + -IA-5/I 2 + Co(a). 

dt -(-4 
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By Gronwall’s inequality, we deduce, for t > to, 



This implies 



t—t oo & IS 


Consequently, 


limsup |u(t )| 2 < — (ca 3 L 3 + — |A 2/| 2 + C 0 (a)) + 2|0 ai£ | 2 



UL 


< —(ca 3 L 3 4—|A 2 f\ 2 + Co(a)) 


a 


u 


by virtue of (12.5ft . where a >8(3. 

Finally, if we choose a = 8(3, then for t sufficiently large, u[t) enters a ball in H 
with the radius 


p 2 := c{(3 2 L 3 + ^r|/H /| 2 + u/3 2 L + (3 3 L 3 + —M h L 3 ). (2.16) 

up u 4 


For simplifying the expression of p, we can assume u < 1 and L,/3 ,7 > 1, then 



This gives p ~ L 5 / 2 , p ~ p ~ /3 5 / 2 , p ~ 7 . The proof of Theorem 12.21 is 


complete. 


□ 


Remark 2.4. By virtue of (12.161) . we see that the dispersion in the KBS equation 

(11. ip tends to amplify the size of the absorbing ball. In particular, if we remove pu xxx 
from the KBS, then p ~ L 3 / 2 (instead of p ~ L 5 / 2 for the KBS). This agrees with the 
estimate in Theorem 2 [37] for the Burgers-Sivashinsky equation u t + uu x = u + u^a,. 

Remark 2.5. By switching t and to in (12.15[) . we obtain, for t < t 0 , 



Assume u[t) : R —* H is a global solution outside the global attractor. By the theory 
of attractors (see, e.g. [63]), we know that |w(t)| is not universal bounded for all 
tell. Therefore, F a:e (t) is not uniformly bounded on R. So, we can choose to such 


that F a ^(to) > ±(ca 3 L 3 + £\A 2 /1 2 + C 0 (a)), and for such t 0 , (I2.17P implies that 


lim inf 


log F a}e (t) ^ a 

|t| -4‘ 
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It follows that 


lim inf 

t —oo 


log |«(t) | 

\t\ 



Since a > 8/3 is an arbitrary large number, we conclude that the growth rate of |w(£)| 
as t —> — oo is faster than exponential, if we assume u(t) is a global solution on M 
outside the global attractor. In the next section, we shall show that in fact, for the 
KBS equation (11,11) with v > 0, there does not exist a global solution outside the 
global attractor. 


2.2. Blow-up of solutions backward in finite time. This subsection is devoted 
to completing the proof of Theorem 11.11 We shall show that, if v > 0, then any 
solution u(t) : [0, oo) —> H of (11.11) outside the global attractor cannot be extended 
to a global solution on M. Our argument follows the approach in [50]. First let us 
prove the following lemma. 


Lemma 2.6. Assume v > 0 and f € V'. Let a > 8f3 and R > ^ +C °^ Q)) ^ _ 

Assume that u(t ) : [£i,£ 2 ] —> H is a solution of lil.ll) such that |u(£i)| < 2 R and 
\u(t)\ > R for t\ < t < £ 2 . Then 


^ 30 

£2 ti — —S 
a 


caL 3 

R 2 ' 


Proof. Recall F a ,e(t) = inf^ g T f T (u(x,t ) — 4> a ,e(x + ff)) 2 dx. It follows that 

F a , e (ti ) < 2(|u(£i)| 2 + \( f a , e \ 2 ) < 8 R 2 + ca 2 L 3 , (2.18) 


where we have used (j2.5[) . 

Notice that the inequality (12.141) states ^ + F a)6 (£) < —f |u| 2 + ^\A~^f\ 2 + C 0 (a). 
Integrating between ti and t 2 and using the fact |w(£)| > R for t\ < t < t 2 , we infer 

F a ,e(t2) — F aje (ti) < — R 2 + (— \A 2 f\ 2 + Co(a))^J (t 2 — ti). (2.19) 

Combining (12.18ft and (I2.19j) yields 

F a ,e(t 2 ) < 8 R 2 + co?L 3 + —R 2 + (—\A 2 ff + C'o(a))^ ( t 2 — £ 1 ). (2.20) 

On the other hand, since F a ^{t) = f T (u(x,t ) — 4> a ,e(x + £,*(t))) 2 dx, then by (12.51) 
we have 

F af (t 2 ) > l|ti(i 2 )| 2 - \4> a .,\ 2 > 1-fi 2 - ca 2 C. 


(2.21) 
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It follows from (I2.20j) and (12.211) that 

X -R 2 - ca 2 L 3 < 8R 2 + ca 2 L 3 + ^-|i? 2 + (^|bH/| 2 + C 0 (a))^ (t 2 - *i). 

This implies 

f R 2 + ca 2 L 3 

t2 ~ U ~ f^ 2 -(f|A-i/| 2 + C 0 (a)) 

^ 30 A 2 + ca 2 L 3 
~ aA 2 
30 | caL 3 
= "a + A 2 ’ 

due to the fact \cxR 2 > 4|A _, 5/| 2 + C 0 (a) from assumption. □ 

Now we can complete the proof of the second part of Theorem 11.11 the blow-up of 
solutions of the KBS equation (II. ip backward in finite time, if v > 0. 


Proof. (Theorem 11.11 (ii)) Assume that there exists a global solution u(t) : M —> H 
which does not belong to the global attractor. By Remark 12.51 |n(t)| grows faster 
than exponentially backward in time. Then, there exists a sequence of times 
with tj + i < tj, such that \u(tj + i)\ < 2■ 7+1 R 0 and \ u (t)\ — 2 J A 0 for tj +i < t < tj, 
j = 0,1,2 ,---, where A 0 > 0 will be chosen later. 

By taking R = 2 3 R 0 in Lemma 12.61 we have, if A^~ l R^aj > ^\A~^ f\ 2 + Co(otj) and 
aj > 8/3, then 


OO OO OO r o 

E , \ > 1 \ — "v CQtj Jj 

(tj - t j+ 1) < 30 — + X/ 

„ Lb n 
3 =0 2 


3=0 


3=0 


4 j R 2 0 


( 2 . 22 ) 


Now we carefully select the values of aj and Rq to force the right-hand side of 
(12.22p to be hnite. We choose aj = 2^ 2 ■ 8/3 A,'/ 4 . Recall Co (a) = cva 2 L + ca 3 L 3 + 
c/3a 2 L 3 + ^-a 5 L 5 . Then, C 0 (a j ) = c(Vv(3 2 rI ,2 L + 2 3 i/ 2 p 3 R 3 0 /4 L 3 + 2^ 3 R l J 2 L 3 + 
1 i 2 5 ^ 2 /3 5 R^ 4 L 5 ). In order to match the requirement 4 3 ~ l Rf a j > ^\A~^ f\ 2 + Co(aj) 
with aj = 2 3 1 2 ■ 8/3R^ 4 , for all j = 0,1, 2, ■ • •, i.e., 




/| 2 • 2~i j + c(2-^ j vPRlL + 2~ 3 P 


-■k n2 


2 rIl 3 + 


-fj 


p 


2 rIl 3 + 


7 


P 


4 4l 5 ) 


for all j = 0,1, 2, • • •, a sufficient condition is 

Ao>^|bl-M 2 + c-max (3 2 L 3 , ^/3 4 L 5 j + L 


(2.23) 
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With oij = 2 J / 2 • 8(3 and Ro satisfying (i2.23j) . we obtain from (12. 22ji that 

<r <t - t \ < 15 w 1 if -a 1 p l 3 

u 1 3 ~ 4 h 2i '^ R T h ^ W 4 A 1 ' 

This can be interpreted as: if |w(t 0 )| > -R 0 , and let the time t goes backwards, then the 
lifespan of the solution backward in time from t — to is shorter than *7 ( ^1/4 + “77?) • 

Therefore, there does not exist a global solution u(t) : M —* H of (II.ip unless u 
belongs to the global attractor. □ 

We remark that the proof of Theorem 11.11 implies that S(t)H is not dense in H 
for any t > 0, where S(t) is the solution semigroup of the KBS equation (II. IK with 
v > 0. This fact is contrast to the Bardos-Tartar conjecture on the 2D NSE [5]. More 
precisely, we have the following result. 

Corollary 2.7. Assume v > 0 and f G V'. For any t' > 0, there exist r(t') > 0 
such that |S'(f)uo| < r(tfor all uq G H, and all t > t', where S(t) is the solution 
semigroup of the KBS equation M.l\) . 

Proof. The result can be directly obtained by the following statement in the proof of 
Theorem 11.11 if u{t) is a solution of (11.11) with |w(£o)| > Rqi then the lifespan of this 
solution backward in time starting at t — to is shorter than + -pyr )• Now, 

given any t' > 0, we choose r(t') sufficiently large so that C ( ^^174 + 7^774) < t'. 
Then, for any u 0 G H , the energy |S'1 < r(f') for all t > f', i.e., the trajectories 
of all solutions stay in a ball of radius r{t') in F[ for all t >t'. □ 


3. Additional Examples 

In this section we discuss the backward behavior of another three dissipative evo¬ 
lution equations, in order to shed light on different types of backward blow-up mech¬ 
anism. We also mention some open problems as well. 

3.1. Damped driven nonlinear Schrodinger equation. I 11 this subsection we 
consider the ID damped driven nonlinear Schrodinger equation on torus (circle) 

iu t + u xx + \u\ 2 u + iXu = /, x G T, (3.1) 

where A > 0, so that iXu is a weak damping. The conservation form of this equation 
(A = 0 and / = 0) has been extensively studied as a fundamental equation in modern 
mathematical physics (see, e.g. 0 129])• For the case A > 0 and / G Lr per ( T), it is 
shown in [3l] that (13.ip is globally well-posed in H per (T) for t G R and possesses a 
weak attractor in H per ( T), and this result has been improved in [67] that the weak 
attractor is in fact a global attractor A in the usual strong topology sense, and that 
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it has finite Hausdorff and fractal dimensions. Also, it is shown in [35] that the 
global attractor A of (13.ip is smooth (i.e. C°°) provided the forcing term is smooth. 
Furthermore, it is proved in [60] . by employing the Gevrey class technique, that A 
is in fact contained in a subclass of the real analytic functions provided the driving 
term / is real analytic. Finally, a recent work [40] shows that A is embedded in 
the long-time dynamics of a determining form (an ODE), and there is a one-to-one 
identification with the trajectories in A and the steady states of the determining form. 

The backward behavior of the solution u[t) of (13.11) . with u 0 G Hf er (T), can be 
described as follows: if u(t) does not belong to the global attractor, then as t —> 
—oo, the L 2 - norm of u{t) grows exponentially fast, while the H 1 -norm grows at most 
exponentially in time. In particular, we have the following: 

Theorem 3.1. Assume f G L 2 er (T). Let u(t ) be a global solution of ( Iff. 1\) in H 3 er (T) 
for all i G 1, with the initial data uq G H 3 er {T), such that u does not belong to the 
global attractor. Then 

c ie - 2 Xt < \u(t)\ < C 3 e~^ xt , as t —>■ — oo, (3.2) 

where c 3 and C\ depends on |u 0 |, |/| and X. In addition, the behavior of the energy 
Eft) backward in time satisfies 

Eft) := \u x ft )\ 2 — - \u(t)|| 4 < C 2 e~ 3Xt , as t —>■ — oo, 

where C 2 depends on ||u 0 ||? |/| and L. Furthermore, 

||w(t)|| < C 3 e~% xt , as t —>• —cxo, 

where C 3 depends on |u 0 |, |/| and A. 

Proof. See the appendix. □ 

Finally, we point out that our purpose of discussing equation (13.11) is to introduce 
the backward behavior of solutions for a typical nonlinear dispersive equation with a 
linear damping: the solution can be extended to all t G R and it grows exponentially 
backward in time. 

3.2. Complex Ginzburg-Landau equation. In this subsection we consider the 
complex Ginzburg-Landau equation: 

u t — (a + bi)u xx — 5u + (a + /3i)\u\ 2 u = 0, in T, (3.3) 

where a > 0, and b,6,a,/3 are all real-valued. Roughly speaking, this is a strongly 
dissipative version of (13.1|) . 

First of all, we notice that, if a = a = 0, then the equation (j3.3[) reduces to a 
nonlinear Schrodinger equation with a linear perturbation, which is globally well- 
posed for t G M. We have discussed this case in section 13.11 

Next, we study a typical parameter regime for equation (I3.3p : a > 0, a > 0. In 
this case, equation (|3.3[) is globally well-posed forward in time and has nontrivial 














BACKWARD BEHAVIOR OF DISSIPATIVE EVOLUTION EQUATIONS 


15 


asymptotic behavior, which has been extensively studied in the literature (see, e.g. 
[2, [231 3D, SB] and references therein). The backward behavior of the solution to (13.3jl 
in this parameter regime can be investigated straightforwardly. Indeed, we have the 
following: 

Theorem 3.2. Assume a > 0 and a > 0. Let u(t ) : [0, oo) — > L 2 {T ) be a solution of 
( OP]) which does not belong to the global attractor. Then u(t ) cannot be extended to 
a global solution of \3.3\) for all ( 6 K, i.e., any extension of u(t) backward in time 
blows up in finite time. 


Proof. We multiply (j3.3j) with u and integrate over T, and then take the real part. It 
follows that 

1 d 


2 dt 


\u\ 2 + a\u x \ — 8\u\ 2 + a|u|]U = 0. 


Since a > 0, we obtain 


1 d 

2 dt 


u\ 2 < S\u\ 2 — a\u\ A L 4 < 5\u\ 


a, 


\u\ 


by Cauchy-Schwarz inequality. By setting s = —t and y = |w| 2 , one has 

d'!J . 2a 2 

Ts >-2 6 V + T y . 

Hence, if y(so ) > —, then y(s ) blows up at some finite time s = Si > s 0 . 


(3.4) 


□ 


Theorem 13.21 states that any solution u to (13.3p off the global attractor blows up 
backward in finite time with respect to the L 2 (T) norm, provided a > 0, a > 0. We 
stress here, since we drop the viscosity in the above estimate, the backward blow-up 
in the case is solely due to the nonlinearity, and it does not result from the presence of 
the viscosity. This is a totally different mechanism, contrast to the backward blow-up 
phenomenon of the KBS equation (11.11) . which is caused by two forces: the cascade of 
energy to small scales by the nonlinearity uu x , as well as the amplification of higher 
modes by the viscosity. In a word, our purpose for discussing equation (j3.3[) is to point 
out a backward blow-up mechanism for certain dissipative PDEs, which is completely 
due to the effect of the nonlinearity. 

Another typical parameter regime for (I3.3[) is that a > 0, a < 0. Under this 
scenario, equation (I3.3j) is locally well-posed forward in time, and by referring to 
the energy identity (13.4j) . we see that the nonlinear term acts as a source forward 
in time, which can lead the solution to blow up in finite time [5] (see also [52] [53J 
for other results concerning finite time blow-up of similar equations). On the other 
hand, for the CGLE defined in the whole space M , the blow-up phenomenon has 
been carefully discussed in [13] [55] [62], etc. However, since a < 0, we see that the 
nonlinearity tends to damp the energy as the time goes backwards, thus the life span 
of solutions backward in time might be a subtle problem, which deserves a future 
study. 
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3.3. Hyperviscous Navier-Stokes equations. We have mentioned in the intro¬ 
duction that, in studying the backward behavior of the 2D periodic NSE [IS], a pair 
of orthogonality relations of the form ((w • V)u, v ) = 0 and ((w • V)w, A u) = 0 is criti¬ 
cal. Also, similar orthogonality relations are valid for the 2D periodic Navier-Stokes-o 
model (65] . However, consider the 2D periodic hyperviscous NSE 

{ u t + z/A 2 u + ( u ■ V)u + Vp = /, x G T 2 , 
div u = 0, (3.5) 

0) = M 0 , 

with mean value zero assumption on Uq, f and u. Notice that the scalar product 
of the nonlinear and linear terms of (13.5p . i.e. ((w • V)w, A 2 w), does not necessarily 
vanish, so the technique presented in [T 8 ] seems not directly applicable to (13.51) . Thus, 
the backward behavior of equation (j3.5j) is an interesting problem to investigate. 

One may first consider a simpler case / = 0. In the absence of external forces, 
the energy satisfies |4|w| 2 + zz|Am | 2 = 0, then by Poincare inequality, we see that 
the energy |w(f)| decays to zero at least exponentially in time, as t —)■ oo, and the 
system has a trivial attractor. In fact, one can show a stronger result: the decay 
forward in time is exactly of exponential type, i.e., \u(t)\ cannot decrease faster than 
exponentially as t —>■ oo. Concerning the 2D or 3D NSE on bounded or periodic 
domain without forcing, the same exponential decay result is also valid, which was 
proved by Foias and Saut in [28], and more precisely they showed that the ratio of 
the enstrophy over the energy converges to an eigenvalue of the Stokes operator as 
t —» oo. A similar result also holds for the solution of (13.51) with / = 0, i.e., the ratio 
|Aw(t)| 2 /|M(t )| 2 is convergent to an eigenvalue of the bi-Laplacian as t —> oo. 

In order to state the following result, we define the space V. :=closure of V in 
Lperi T 2 ) where V = {(p = trigonometric polynomials with values in K 2 : div <p = 
0, f T 2 pdx = 0}. The well-posedness of (13.51) for t e [0, oo) with u 0 G T-L is classical. 

Theorem 3.3. Let / = 0 and uq G "H, then the solution u of \3.5\l satisfies 

e —&(t—t 0 )| w (£ o )| < l'u(t)! < e _ 1 ' A i ( - t-to ^|w(to)|, for t > to > 0, (3.6) 

where b = v exp jcz^A^mto) | 2 | an d Ai = (2n/L) 2 . 

Proof. The proof follows the approach in [28]. For the sake of completion, we provide 
the proof in the appendix. □ 

Remark 3.4. Assume that there exists a global solution u(t), ( 6 1, of (13. 5 p with 
/ = 0. By switch t and t 0 in (13.61) . we obtain 

\ u (t)\ > e uX ^ to ~^\u{t 0 )\, for — oo < t < t 0 < oo. 

This shows that the energy \u(t)\ of a nonzero global solution u grows at least expo¬ 
nentially fast as the time t goes backwards. Thus, there are three possibilities of the 
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backward-in-time behaviors of the solution u, that is, |u(i)| increases exactly of ex¬ 
ponential type as t —» — oo, or faster than exponential, or blow up backward in finite 
time. Here, concerning the 2D periodic hyperviscous NSE (13.5p . we state a similar 
conjecture to the one of Bardos-Tartar [5j: S{t)l-L is dense in % for every t > 0, where 
S(t) is the solution semigroup of (13.5p . In fact, we can propose a slightly stronger 
conjecture for (j3. 5p : the set of all u$ G "H for which S(t)uo is a global solution for 
t G R is dense in %. 

Finally, we point out that, for the 2D periodic hypervicous NSE (13.5p without 
forcing, i.e., / = 0, there exists a class of initial data u 0 for which the solution S(t)u 0 
of (13.5p is global, for all t G K, such that the energy grows exponentially fast 

as t —* — oo. In order to construct these explicit solutions we write (13. 5p with / = 0 
in the equivalent vorticity-stream formulation (see, e.g. 0) 


Lu t + z/A 2 cu + J(0, A0) = 0, x G T 2 , (3.7) 

where the vorticity cu := curl u, 0 is the stream function satisfying u = := 

(—'0a; 2 , Vhi ) tr , and the Jacobian J(0, A0) := ^ x , 2 . Let us define, for any 

fixed p G N, 


ipo(x) 


E 

keZ 2 , |k| 2 =/3 


r ,27rk \ , . /27rk 

[a k cos(—— ■ x) + b k sm( — — • x)J, 


x G T 2 , 


where a k and b k are arbitrary real numbers. Notice that 0 O is a steady state of the 
2D Euler equation uj t + J(0, A0) = 0 on T 2 , due to the fact J(0 0) A0 O ) = 0 (see, e.g. 
H). Inspired by these steady state solutions of the Euler equation, we see that 

■0OM) = e~ 16n4 P 2ut/L4 Mx) 

is an explicit solution of (13.71) . since for this particular stream function, one has 
J(0, A-0) = 0, and the vorticity 

u(x,t) = A0(x, t) = e _ 167f4 ^ 2 Ut/Ij4 0o (x). 

Also, we can calculate the velocity 

.. _ vvLy, _ ( Sk=(fc 1 ,fc 2 )eZ 2 , |k| 2 =^[ fl k^2 sin(^ ■ x) - b k k 2 cos(^A . x )] \ 
VEk=(ifc 1 ,fe2)ez 2 ,|k| 2 =^[-«kfc 1 sin(^-x) + 6 k fc 1 cos(^-x)] > / ' 

Since (13.7p is equivalent to (13. 5 p with / = 0, we have that 

u(x,t) = V x 0 = e-^^uo 

is an explicit global solution of (13. 5 p with / = 0, for all f el. Obviously, the L 2 -norm 
\u(t)\ = e - 16n4 P 2ut / L4 I^Q| grows exponentially fast backward in time. In other words, 
there exists a family of initial data uq, constructed above, for which the solution 
S(t)uo of (13.5p is global, for all f el, with the L 2 -norm growing exponentially fast 
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backward in time. However, it is an open problem whether that the set of initial data, 
for which the solution of (13.5j) can be extended to all t G M, is dense in "H. 


4. Discussion 

In this section we would like to provide some discussion for the purpose of un¬ 
derstanding the different backward behaviors of the 2D periodic NSE and the KBS 
equation (11.11) . In section I4~T1 we will consider a special case of the KBS equation: the 
viscous KdV equation, and demonstrate that a cnoidal-like wave solution (periodic 
soliton) blows up backward in time. In section 14.21 we employ the energy spectra 
and argue, using Kolmogorov turbulence theory, in an attempt to explain why the 
solutions of the 2D periodic NSE can be extended to a global solution, for all t 6 1, 
for a rich set of initial data, but any solution outside the global attractor of the 
KBS equation (11.11) . with u > 0, must blow up backward in finite time. Our discus¬ 
sion is intended to shed more light on the nature of this phenomenon using physical 
non-rigorous arguments. 

4.1. Physical interpretation of backward behaviors of perturbed KdV. We 

aim to discuss the backward blow-up of the KBS equation ( 11 . 11 ) using perturbation 
arguments on explicit solutions. For the sake of simplicity, we set f3 = 0, / = 0 and 
v > 0 in (HU), the equation becomes a viscous perturbation of the KdV equation 

Ut 4“ mix T 'y'Uxxx 0, in If — [ L/2, Tjj 2], (4.1) 

with the spatial mean value zero assumption on u. We have shown that all of the 
nonzero solutions of (14.11) blow up backward in finite time. It is well-known that the 
existence of the soliton solutions for KdV represents a balance between the dispersion 
and nonlinearity. Adding a viscous term to the KdV equation affects such balance, 
the viscosity dissipates the kinetic energy and tends to diminish the soliton forward 
in time. However, backward in time, the viscosity acts as a source that may intensify 
the soliton, and cause the deformation of its shape, leading to higher amplitude and 
shrinking width. Thus the blow-up of the soliton backward in time may result from 
the amplitude of the perturbed soliton approaching infinity. 

To demonstrate this argument more clearly, we shall look at an explicit periodic 
soliton solution of the KdV equation and employ a standard perturbation argument. 
Let us consider the classical KdV equation 

u t + uu x + u xxx = 0, in T. (4.2) 

It was discovered by Gardner, Greene, Kruskal, and Miura in their seminal work 
[30) (31] that, if u is a solution of the KdV (14.21) . then the eigenvalues of the Schrodinger 
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operator d xx + are independent of time. Specifically, suppose that there exist A (t) 
and a function 0(x,f) G H 2 { T), with f T ift 2 (x)dx = 1, satisfying 

{dxx + -u)^ = X'i/j, in T, (4.3) 

6 

then A(f) is a constant in time. This can be seen readily by a formal argument [30J. 
Indeed, a straightforward manipulation of (14.31) implies 

d 1 

— ( 'IpxR ^Rx) 'Ip ri U t r U j xxx)') i ( 4 . 4 ) 

OX o 

where R(x,t) := 0 t ~ \ u x^ + (|m + 4 A)02,. Since u is a solution of (14.2p . then (j4.4j) 
is reduced to 

— i/jR x ) — —'0 2 At- (4.5) 

Integrating (14.5[) in T, by using the periodicity of -0 and w, as well as the square 
integrability of -0, we infer that X t = 0, i.e., A(f) is invariant in time. 

By employing the technique mentioned above, a perturbation of the KdV equation 
can be analyzed. The method is standard and can be found in [U 124) 02] . Indeed, let 
us consider a perturbed KdV: 

u t + uu x + u xxx = eg, in T, (4.6) 

where e > 0 is small and q{x, t ) may depend on u. The goal is to find the evolution 
of the eigenvalues of the operator d xx + provided u is a solution of (14.6lh To this 


end, we consider the equation 



{dxx + \u)‘$ = A0, 

6 

in T, with / 0 2 (x)dx = 1. 

J T 

(4.7) 

Exactly the same as (|4.4|, we have 



d ~ ~ — 

-T^i'lpxR 'Ip R'x') ~ 

~lp (\ t — {lif “I - r U>'U' X “ 1 “ 'U'xxx)') 5 

0 

(4.8) 


where we define R(x, t) = ijj t — | - u x '0 + (|u + dA)^. By substituting (14. 6 p into (14.8p . 
it follows that 

p-{$ x R - $R X ) = 0 2 (A t - \eq). (4.9) 

OX 0 

Unlike the situation of the classical KdV equation (14.2p . here an eigenvalue A(f) is no 
longer invariant in time. Indeed, by integrating (14.9p over T and using the periodicity 
of 0 and u, it follows that 

A t f 0 2 dx = -6 [ i[) 2 qdx. (4-10) 

J T 6 Jf 
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Since J T ip 2 dx = 1, then 

A t = \t [ ftqdx. (4.11) 

6 J T 

Note that (14.lip shows that an eigenvalue A (t) may evolve slowly within a short time 
provided e is small. 

Next we look at a typical explicit periodic solution of the KdV equation defined in 
the torus T = [--§,§]• Take the initial data to be uq = 12m 0 /o cn 2 (/ 0 a;, m 0 ), where 
mo G (0,1) is the elliptic parameter for the Jacobi elliptic function cn(;r,mo). Here, 
we require mo and lo to satisfy the relation 

1 f 2n 1 

*o = 7 / 7 . (4.12) 

L J o y/1 — mo sin' 0 

so that the period of uq(x) is L , due to the definition cn(x,m 0 ) = cos(0) where 
x = f n ' ■. — 1 --~ =dd. With this initial data u 0 , the unique solution of the KdV 

u y 1-mo sin 2 0 

equation (14. 2 p is 

u(x, t) = 12m 0 /o cn 2 (Z 0 ((£ — c 0 t),m 0 ), (4-13) 

where Co = 4(2mo — 1)Zq (see, e -g- ; [211 51 j). This is the so-called “cnoidal” wave 
solution (periodic soliton) derived by Korteweg and de Vries in their 1895 paper |45] 
in which they proposed the KdV equation. Concerning the spectral problem (14.31) 
with the cnoidal wave solution u given in ( 14 .1 3h . we find an eigenvalue 

A = (2m 0 - l)/ 0 2 , (4.14) 

with the corresponding eigenfunction 


V’OM) = 


a 


cn (Z 0 (a;-c 0 t),m 0 ), 


(4.15) 


mo 


where C mo = r^'o/^ cn 2 (a;, mo )dx, where K mo = f^ n 1 =d6 is the period of 

the elliptic function cn(x,mo). With this definition of ip, one has J T ip 2 (x)dx = 1. 

For the same initial data uo = 12m 0 /o cn 2 (/ 0 a:, mo), we now investigate the solution 
u of the perturbed KdV equation (14.61) provided u(x, 0) = u 0 - Since e > 0 is very 
small, we expect the solution u of (14.6p to be a small perturbation of the KdV solution 
(14. 130 . at least for a short time, in the sense that, 


u(x, t ) ~ 12ml 2 cn 2 (l(x — ct),m), 


(4.16) 


where l(t) = \ n , 1 ==d9 and c(t) = 4(2 m(t) — l)l 2 (t), with 0 < m(t) < 1, 

L JU yj l—m{t) sin 2 6 

such that 1(0) = lo and m(0) = mo- As a result, for the spectral problem (14.71) . it is 
expected that there exists an eigenvalue A with the corresponding eigenfunction ip, 
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which are small perturbation of A and i/i, defined in (14.141 and (14.151) . respectively, 
he., 


A (t) r; (2 m(t) — 


(4.17) 


with the corresponding eigenfunction 


4>{x,t) R 





(4.18) 


where C m (t ) = 


J—K ^/2 cn2 ( x , m (t))dx, where K m (t ) 



_l_ 

^/l—m(t) sin 2 0 


dO is the 


pe¬ 


riod of the elliptic function cn (x,m(t)). 

In order to find out how the shape of the cnoidal-like wave solution u changes in 
time for the perturbed KdV equation (j4.6[) with small perturbation eg, we study some 
typical perturbations. 


4.1.1. The damped KdV equation. First set q = — u in (14.61) . i.e., we consider the KdV 
equation with a small linear damping 

u t + uu x + u xxx = -eh, in T. (4.19) 


Equation (j4.19j) is a special case of a more general equation u t + uu x + u xxx = — T{t)u , 
which models a wave moving over an uneven bottom .'19. 0D Db By (14.161) , q = 
—u r — 12ml 2 cn 2 (Z(a; — ct),m), then we substitute this along with (I4.17P and (14.181) 
into (14.lip to obtain 


—2 eml 2 r K ™/ 2 

d t [(2m — l)l 2 } r—— - / cn 4 (x, m)dx. (4.20) 

V m J — Km/2 


In order to see the evolution of l(t) more clearly, we consider the situation in which 
the elliptic parameter m is very close to 1, so that m is almost invariant in short 
time compared to the change of /, and then the Jacobi elliptic function cn(x,m) can 
be approximated by the hyperbolic function sech(a;). Under such scenario, we obtain 
from (I4.20p that 


d t l ~ —eCl, with 1(0) = l 0 , 


(4.21) 


where lo is given in (14.12p . and C = (J R sech 4 (a;)(ia;) (J K sech 2 (a;)(ix) 1 = |. Notice 

that (I4.2ip implies l(t) Ri e~ eCt lo , i.e., forward in time, l(t ) decreases to zero exponen¬ 
tially fast. On the other hand, backward in time, l(t ) increases with an approximately 
exponential rate, i.e., the amplitude 12ml 2 of the cnoidal wave (14.161) grows exponen¬ 
tially with an exponentially fast shrinking width ~ 1 //. Moreover, the propagation 
speed c(t) = 4(2 m(t) — l)/ 2 (t) of the cnoidal wave also increases exponentially back¬ 
ward in time. Furthermore, by (14.16P it follows that the energy \u(t)\ 2 ~ l 3 (t), which 
increases exponentially fast backward in time. Notice that this finding is consistent 
with the L 2 energy estimate of (14.19p . i.e., \u(t)\ 2 = e~ 2et \uo\ 2 . In addition, from 
(I4.16P we see that \u x (t)\ 2 ~ l 5 (t ) also tends to grow exponentially backward in time. 
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4.1.2. The viscous KdV equation. Next we set q = u xx in f!4.6j) . i.e., we consider the 
viscous KdV equation 

u t + uu x + 

^XXX eu xx , in T. (4.22) 

By using (I4.16|) . we calculate 

q = u xx ~ —24m/ 4 [3mcn 4 (/(x — ct),m) + (2 — 4 m) cn 2 (/(x — ct),m ) + m — 1]. 
Substituting this formula along with (14.171) and (14.181) into (14. 11 j) gives 

4 rm / 4 t'Km/ 2 

d t [(2m — 1 )l 2 } ~-—— / cn 2 (x,m)[3mcn 4 (i,m) 

C m J —Km/2 

+ (2 — 4m) cn 2 (x, m) + m — l]dx. (4.23) 

Analogously as above, we consider the case that the elliptic parameter m is very close 
to 1 so that m is almost invariant in short time compared to the change of l, and 
then cn(x, m) can be approximated by sech(x). Thus, (14.231) implies 

d t l « —eCil 3 (4.24) 

where C\ = 2 (J R sech 4 (x)[3sech 2 (x) — 2 }dx) (J R sech 2 {x)dx) 1 = pr. It follows that 
/(t) decreases to zero forward in time. In addition, we see from (j4. 24[) that l[t ) tends to 
approach inhnity backward in finite time, i.e., the amplitude 12ml 2 of the cnoidal wave 
(I4.16j) grows very fast with a rapidly shrinking width ~ l/l and a fast accelerating 
wave speed c = 4(2m — 1)Z 2 , which may lead to a backward blow-up. Also, we obtain 
from (I4.16j) that the energy \u(t)\ 2 ~ l 3 (t), tending to blow up backward in time. 
This finding indicates the mechanism of singularity formulation, and it is consistent 
with the rigorous mathematical result in Theorem 1 (ii) that the energy of the viscous 
KdV (14 .22 p blows up backward in time. 

4.1.3. The viscous Benjamin-Bona-Mahony type of equation. Finally, in order to test 
the robustness of the asymptotic argument presented in this section, we would like to 
study another perturbation of the KdV equation. In particular, we set q = u xx + u xxt 
in (14. 6p . and then the equation becomes 

Ut T uu x ~\~ u xxx V U'xxt)i hi T. (4.25) 

Equation (j4.25[) is a viscous Benjamin-Bona-Mahony (BBM) type of equation. The 
original BBM equation, u t + u x + uu x — u xxt = 0, was proposed as a modification of 
the KdV equation for modeling long surface gravity waves [4]. The energy identity of 
(14.251) reads 

-_ r(|h | 2 + e|Aj;| 2 ) = —e|it x | 2 . 

2 di 1 1 1 1 11 

This implies, with the Poincare inequality, that the energy |u| 2 + e|'U 3; | 2 decreases to 
zero exponentially fast forward in time, while grows exponentially fast backward in 
time. Now we use our perturbation argument to study the evolution of the shape 
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of the cnoidal wave for (14.25j) . Similar to arguments introduced above, we assume 
that the elliptic parameter m is very close to 1 so that m is almost invariant in short 
time compared to the change of /, and in this case, cn (x,m) can be approximated by 
sech(x). Then, after straightforward calculations, we obtain 

-2 l t « 2 eC 2 l 3 + 32eC 3 l 6 + l t (8eC 2 l 2 - 8 eC 4 l 2 + 64 eC 3 lh), 

where 

C 2 
C 3 
c 4 


= / sech 4 (x)[3sech 2 (a;) — 2 ]dx = — ; 


15 


= / sech 4 (x)tanh(a:)[3 sech 2 (a:) — 1 ]dx = 0 ; 


= / sech 4 (a;)tanh(a;)[3sech 2 (a:) — l]xdx = 


It follows that 


d t l 


8 el 3 

20 d 2 + 15' 


On one hand, this implies that l(t) decreases to zero forward in time. On the other 
hand, backward in time, for large l, the amplitude 12mZ 2 of the cnoidal wave (14.161) 
grows approximately exponentially fast with exponentially fast shrinking width i A 
which is consistent with exponential growth of the energy backward in time. 


4.2. Energy spectra. In this section, we would like to argue that the various back¬ 
ward behaviors of different dissipative equations are connected to their energy spectra 
as well as the Kolmogorov turbulence theory. We stress that our arguments are mainly 
physically oriented, rather than a rigorous mathematical treatment. 


4.2.1. The viscous KdV equation. In the previous discussion, we have given some 
evidence that, if we take the cnoidal wave to be the initial value of the viscous KdV 
equation (14.11) . then as the time goes backwards, the amplitude ~ l 2 of the cnoidal 
wave grows fast with its width ~ 1 shrinking rapidly, and as a result, the solution 
may get close to a Dirac delta function. Since the Fourier transform of a Dirac delta 
function is a constant, we realize that the energy spectrum El D of the solution u, 
defined as El D = \uk\ 2 , for k <E Z\{0}, is approximately invariant with respect to k, 
for the spatial scales j|t sufficiently larger than the width ~ j of the cnoidal wave (see 
(S3]). Now, by setting s = —t, we consider the energy identity of the viscous KdV 
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(14.lj) . backward in time, 

= y\ u x{s )\ 2 = u k 2 E l D ( s ) 

fcez\{ 0 } 

«■'(c E fc2 + E *? E l D ( s n ■ ( 4 - 26 ) 

y o<|fc|<fc max (s) \k\>k max(s') J 

Notice that k max (s) ~ l(s)L, where L is the length of the domain T. Recall we have 
argued in subsection 14.1.21 that l(s) tends to blow up in finite time, thus k max (s ) 
increases very fast, which leads to a fast growth of -j-\u(s)\ 2 due to (14.261) . ending up 
with an energy blow-up in finite time. 

4.2.2. The Burgers equation. The KBS equation (II. ip is also related to the Burgers 
equation. In fact, by setting (3 = 7 = 0 and / = 0 in (II. ip . it reduces to the Burgers 
equation u t — vu xx + uu x = 0, and we have proved that all of its nonzero solutions 
blow up backward in finite time. It is well-known that the Burgers spectrum behaves 
like El D ~ k~ 2 (see my An energy formula analogous to (j4.26j) indicates the energy 
tends to grow very fast as the time goes backwards, leading to a finite time blow-up. 

4.2.3. The 2D Navier-Stokes equations. Next we attempt to discuss the relation be¬ 
tween the Kolmogorov turbulence theory and the backward behavior of 2D NSE. The 
Richardson energy cascade theory [63j states that the energy transfers among eddies 
with similar sizes only, and the rate of the energy injection at large scales is equal on 
average to the energy dissipation rate at small scales, so that a statistically steady 
turbulent state forms. Far away from the source and sink, there is an inertial range, 
in which the turbulence properties only depend on the energy cascade rate. Based on 
a simple dimensional argument, in 1941, Kolmogorov and Obukhov [ 431 [441, 59] de- 

2 5 

rived a celebrated result of the energy spectrum for the 3D turbulence ~ C e e 3/j“3 
provided the wavenumber k > 0 is in an inertial range, and e is the energy cascade 
rate (equal to the energy dissipation rate). 

On the other hand, the 2D incompressible ideal flow conserves two quadratic quan¬ 
tities, energy and enstrophy, and thus the 2D turbulence possesses a dual cascade 
behavior [46]: inverse energy cascade and direct enstrophy cascade, i.e., the energy 
is transferred from small to large vortices, while the enstrophy is transferred to small 
scales. By a dimensional argument, the inverse energy cascade spectrum is identical 
to the 3D turbulence spectrum mentioned above, but the enstrophy cascade spec¬ 
trum reads E K ~ C^^k - 3 , which is called the Kraichnan spectrum [46], where r/ is 
the enstrophy cascade rate, and the isotropic spectrum for each n > 0 is defined by 
E k := 27w|m k | 2 , with the isotropic assumption \u K \ = \u^\ for all set 2 with |k| = k. 
A more precise version of the Kraichnan spectrum with a log correction is suggested 
in m- The Kraichnan dual-cascade picture was recently confirmed numerically in 
[8]. For the turbulent flow in the periodic domain, the energy spectrum E K can also be 
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defined using the following way (see m : it is assumed that there exists E K such that 
the total energy between the wavenumbers k' and k", i.e. J^kez 2 fc'<|k|<fc" l^k| 2 , can 

rk" 

be approximated by the integral J fc , E K dn, for wavenumbers k < k" in the inertial 
range. 

Now, assume the energy is injected near wavenumber kf and dissipated at very small 
wavenumbers kf and at very large wavenumbers k + kf, and there are neither 

forcing nor dissipation at wavenumbers n such that < hi < kf or kf < hi < k + , 
which are the inverse and the direct cascade inertial ranges, respectively. Under such 
scenario, if we set s = — t, then the energy identity of the 2D NSE, backward in time, 
reads 




V 


H 2 E K dhi 



fk+ 


n z E K dn + C e 63 I KadK + C^rjs I k l dn + 


' k— 



(4.27) 


Since in general the lower modes are not a main source for the energy blow-up, we 
focus on the higher modes hi > kf. By a simple dimensional argument, the enstro- 
phy cascade inertial range is expected to extend up to k + rs-/ (r)/v 3 ) s [46]. As a 
result, the third integral in (14.27ft can be estimated as follows: C v 7 /§ K~ l dn ~ 

C r? ?u(| ln(?//z/ 3 ) — In (kf)), which is controlled by the enstrophy cascade rate. Con¬ 
cerning the last term in (I4.27|) for larger wavenumbers n > k + , we notice that, as 
the time goes backwards, the viscosity amplifies the enstrophy on high wavenumbers, 
thus the enstrophy tends to cascade inversely backward in time, moving from small 
to large scales, so that the enstrophy on the small scales might not grow too rapidly. 
Consequently, the energy increasing rate -^\u(s)\ 2 is possible to be relatively slow, 
and the finite-time backward blow-up of 2D NSE may be prevented, for initial data 
with energy concentrated on low modes. 

Our argument above is consistent with the result in [18] , which states that there is a 
rich set of initial data in the energy space % for which the solution of 2D NSE can be 
extended to a global solution for all tGl. More precisely, it is shown in PS] that, if 
Po G P n H, for some n , then there exists a global solution S(t)uo such that P n Uo = Po, 
with \Q n u 0 \ < max{2|/|/z/Ai,7n /2 |p 0 |} where = (A n+i + A n )/(A n+ i - A n ), and 


lim sup 

£—>•—00 


\S(t)uo\ 2 


< 


An + A n +1 

2 


where P n is the projection onto the lower modes \k\ < n, Q n = I — P n , and A n is 
the nth eigenvalue of the linear Stokes operator in 2D periodic NSE. This indicates 
that, if |Q n iio| is controlled by \P h uq\ in an appropriate manner, then the ratio of the 
enstrophy to the energy can be bounded as t — oo, which may be a result of the 
inverse enstrophy cascade and the direct energy cascade, backward in time, from the 
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Kraichnan dual cascade picture for 2D incompressible flow. Also, it is easy to see 
from the energy identity that, if the ratio of the enstrophy to the energy is uniformly 
bounded for all negative time, then the solution is global for all t £ M, growing 
exponentially as t —* — oo, provided it does not belong to the global attractor. 

In sum, the energy and enstrophy cascade as well as the Kolmogorov spectra may 
reveal the underlying mechanism for the different backward behavior of the KBS 
equation (11.11) and the 2D NSE. In particular, we would like to emphasize that there is 
a huge difference between the KdV spectrum E\ D ~ constant and the 2D turbulence 
spectrum E K ~ n~ 3 . For more material on Kolmogorov turbulence theory, please 
refer to monographs ®mmm- 


5. Appendix 

We prove Theorem 13.11 and Theorem 13.31 in the appendix. 

5.1. The proof of Theorem 13.11 

Proof. The proof consists of two parts: the first part is devoted to the L 2 -estimate, 
and the second part is devoted to estimating the energy as well as the i7 1 -norm. 

(1) L 2 -estimate 

Multiply (13.1|) by u and integrate on T, 

i [ u t udx - \u x \ 2 + |u|£ 4 + i\\u\ 2 = [ fudx. (5.1) 

J T J T 

By taking the imaginary part of (15. 1|) . we obtain 

\j t \ u \ 2 + A ^ 2 = Jj udx - ( 5 - 2 ) 

where we used Cauchy-Schwarz inequality and Young’s inequality. 

It follows that 

-||m| 2 + 3A|m| 2 > -1|/| 2 . 

By using Gronwall’s inequality, we have 

\u{t)\ 2 < e" 3At |u 0 | 2 + ^l/l 2 i e ~ 3M ~ X ) > for t < 0. (5.3) 

On the other hand, (15.21) also shows 

t±^ + \\u? = l™fjud X <\\u? + ±\f\\ 


>UAM 2 <i|/| 2 . 


Hence, 
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Again, by Gronwall’s inequality, we obtain 

\u(t)\ 2 >e~ xt \u 0 \ 2 + ^\f\ 2 {l-e~ xt ), for t < 0. (5.4) 

By (15.3p and (15.41) . we infer, if u does not belong to the global attractor, then 

c\e~^ xt < \u(t)\ < Cie~^ xt , as t -* —oo, (5.5) 

where c\ and C\ depends on |u 0 |, |/|, and A. That is to say, the L 2 -norm of u(t), which 
is outside the global attractor, grows exponentially fast as the time t goes backwards 
to negative infinity. 

(2) Estimate of the energy and the H 1 -norm 

Multiplying (13.ip by u t followed by integrating on T, we obtain 


i\u t f 


u x u xt dx + / u 2 uu t dx + iX I uu t dx = I fu t dx. 


Taking the real part gives 


d 

dt 


u x \ — -\u\ L i ) + 2A ( Im I uu t dx ) = — 2 ( Re j fu t dx 


(5.6) 


Now we need to estimate Im J T uu t dx. Indeed, by taking the real part of (15.ip . we 
infer 


Im / uu t dx = —Im / u t udx = \u x \ 2 — \u\ A L 4 + Re / fudx. 


(5.7) 


Substituting (15.7p into (15. 6p . we obtain 


d 


— u 


dt 


-\u\ L 4 + 2A I \u x \ — \u\i4 + Re / fudx = — 2 I Re / fu t dx . 


) 




This can be written as 


—0 + 2A0 = A|w|^ 4 + 2A ( Re j fudx ) , 


(5.8) 


where 0(f) = \u x \ 2 — + 2 (Re f T fudx). 

Note that, \u\ 2 < \u\ 2 l 4 L^ < |rt|| 4 +j, by virtue of Holder’s and Young’s inequalities. 
Thus 

Re fudx > ~\ u \ 2 - i |/| 2 > ^ - \\f\ 2 - ( 5 - 9 ) 

Applying (15.9p on (15.81) . one has 


±t + 2 ^>-xa +w 
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By Gronwall’s inequality, we deduce 

m < e-™m + \{^ + l/l 2 ) (e- 2A ‘ - 1), for t < 0. (5.10) 

Define the energy E(t) := \u x \ 2 — ||u|| 4 . Then 

4>{t) = E(t) + 2 ^Re J fudx^j > E[t) — |/| 2 — \u(t)\ 2 . 

Therefore, by (15. 3ft and (15.101) . one has, for t < 0, 

E{t)< j>(t) + \u(t)\ 2 + |/| 2 

< e- 2J V(0) + 1 + |/| 2 ) (o' 221 - 1 ) + <r"k| 2 + fy\f\ 2 (e~ 3M - 1) + l/| 2 - 

This shows 

E(t) < C 2 e _3/U , as t —)■ — oo, (5.11) 

where C 2 depends on ||w 0 ||, |/| and L. 

Finally, by the ID Agmon’s inequality as well as the Young’s inequality, we deduce 
E(t) = \u x \ 2 — -\u\\i > | u x \ 2 — c|u| 3 ||u|| 

> Ilia,| 2 — - ||u|| 2 — c|w| 6 

> -||u|| 2 — |u| 2 — c|u| 6 . 

Thus 


||w|| 2 < 2 E(t) + 2|u| 2 + c|u| 6 . 
It follows from (15. 5 p and (15.lip that 

_9u 

||w|| < C 3 e 2 , as £ —>■ —oo, 
where C 3 depends on |m 0 |, |/| and A. 


□ 


5.2. The proof of Theorem 13.31 

Proof. The proof follows the idea in [28] . Let / = 0. Recall the space PL is the 
space of all periodic, divergence free, L 2 -functions on T 2 with vanishing mean values. 
Let P be the orthogonal projection on PL in L 2 er (T 2 ) 2 and set Au = —PAu and 
B(v,w ) = P[{y ■ V)w/]. Then the equation (13.5p can be written in an equivalent form 

(5.12) 


u t + vA 2 u + B(u, u) = 0, in T 2 . 
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Since (B(u,v),v) = 0 (see, e.g. [E]), taking the scalar product of (15. 12ji with u gives 


+ = o. 


(5.13) 


Let Ai = (27 t/L) 2 be the first eigenvalue of A = —PA. Then by Poincare inequality 
\Au\ > Ai|u|, we obtain 

+ pA 2 |m| 2 < 0. 

Zj cll 

It follows that 

\u(t)\ A e~ uXl( ' t ~ t °' > \u(to)\, for t > to > 0. (5-14) 

Next we show that \u(t)\ > e _6(t_to ^|u(to)|, for t > t 0 > 0, where b > 0 will be 
specified later. To this end, we set 


Q(t) = 


I MW 


W(t )\ 2 ' 

By differentiating both side of q(t)\u{t)\ 2 = \ Au(t)\ 2 , we obtain 


(5.15) 


d 


—q \u\ + 2 q ( u t , u) = 2 — Au , Au 


It follows that 


1 d 

2 dt 


Q = 


\u\ 


dt 


;(u t ,A 2 u - qu ) 


d 


dt 


\uy 

1 


\u\ 


-(i/A 2 u + B(u, u ), A 2 u — qu) 

-z/| A 2 u — qu\ 2 — y^r(B(u, u ), A 2 u — qu) — u(qu , A 2 u — qu). 


Notice, by (j5.15[) . we have 


. , 9 , \Au\ 2 , .ox lAul 4 , 

(qu, A u — qu) = ——— (u, A u) -j—j— (u, u) = 0. 


\u\ 


\u\ 


Combining the above two estimate and by letting v = A, we have 


1 d 

2 Jt q ~ 


l«l : 

ju\A 2 u — qu\ 2 — - — pr(B(u, u), A 2 u — qu) 


\u\ 


\u\ 


= —u\(A 2 — q)v\ 2 — \u\(B(v,v), (A 2 — q)v). 
Consequently, by Cauchy-Schwarz inequality and Young’s inequality, we infer 


1 d 

2 dt 


q + v\(A 2 — q)v | 2 = —\u\(B(v, v ), (A 2 — q)v) 


v 


<-\(A 2 -q)v\ 2 + -\u\ 2 \B(v,v)\ 2 . 
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It follows that 


4? < -H 2 \ B (v,v)\ 2 = -^\B(u,u)\ 2 . 
at v u\u\ z 


Notice that \B(u,u)\ < \u\ l a\\ 7u\ L ± < c|zt |2 ||rt|| 2 , by virtue of 2D Ladyzhen¬ 
skaya’s inequality \u\n < c|«| 2 ||w|| 2 . Then 

\Au\ 2 ' 


d c 

T + q -~T\ 
dt u\u\ 


|u||| u|| 2 |t4u| < cu 1 A 1 2 |u|||u 


\u\ 


= cv 2 A X 2 |u|||u||g, 


where Poincare inequality has been used. This gives 


q(t) < q(to) exp <! cv 1 A 1 2 / |u| ||u||ds/>, for t > f 0 > 0 


'to 


Now we estimate f* |w|||iz||ds. Indeed, by fl5.13lh one has 


'to 


\Au\ 2 ds = — (|u(t 0 )| 2 - W(t)\ 2 ) . 


(5.16) 


(5.17) 


Therefore, by Poincare inequality, Cauchy-Schwarz inequality, as well as the estimates 
(I5.14p and (I5.17I) . we deduce 


|u|||u||ds < A 2 / |u||Au|ds 


'to 


'*0 


< A 2 


t 

2 


\u\ ds 


to 



t \ 2 

2 . 


\Au\ 2 ds 


to 




<(2v) U, 2 \u(t 0 )f. 

Substituting this estimate into (15.1611 gives 


q(t) < q(t 0 ) exp < cv 2 A x ^|w(t 0 )| 


That is 


I Au(t)\ 2 < 


I Au(t 0 )\‘ 


H*o)| 2 

and along with the energy identity, we obtain 


u (t)\ 2 exp < cv 2 A 1 2 |w(t 0 )| 2 >, 


d 'u(t) P + 2J A “ ((o)I 


dt 


|u(f)| 2 exp jcz/ 2 \ 1 2 |«(t 0 )| 2 | > 0. 


It follows that 


,«(t«,)| a 

u(t)\ > \u(t 0 )\e~ b{t ~ to \ for t>t 0 >0, 
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where b = z/^^exp jcz/ 2 2 |w(t 0 )| 2 j . □ 
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